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Abstract 



We look for differential equations satisfied by the generalized Jacobi polynomials 
{P"'^'^''^ i^)} n-o '^^^'^^ ^-^^ orthogonal on the interval [—1,1] with respect to the weight 
function 

2^.mr(.^^"lr^ + i) ^^ - + + ^^^(^ + ^) + ^^(^ - 

where a > -1, /3 > -1, M > and > 0. 

In order to find explicit formulas for the coefficients of these differential equations we 
have to solve systems of equations of the form 



J2M^)D'Pi"'^H^)=Fn{x), n = l,2,3,. 



1=1 

where the coefficients {Ai{x)}'i^^ are independent of n. This system of equations has a 
unique solution given by 

This is a consequence of the Jacobi inversion formula 
a + P + 2k + l 



f^. (a + /3 + /c + j + l),_j+i 

K—J 

which is proved in this report. 
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1 Introduction 



In ijT^ T.H. Koornwinder introduced the generalized Jacobi polynomials 
which are orthogonal on the interval [—1,1] with respect to the weight function 

r(« + /3 + 2) , 



oo 





2a+/3+ir(a + l)r(/3 + 1) 



1 - x)"(l + x)f + M5{x + 1) + N6{x - 1), 



where a>— 1,/3>— 1, M>0 and > 0. As a limit case he also found the generalized 
Laguerre polynomials \ L'^'^\x) \ which are orthogonal on the interval [0, oo) with respect 

I J n=0 

to the weight function 

—^x^e-- + M5{x), 
r(a + 1) 

where a > — 1 and M > 0. These generalized Jacobi polynomials and generalized Laguerre 
polynomials are related by the limit 

L^'*^(x)=^limP:'/^AM(^l_|^. 

In Q we proved that for M > the generalized Laguerre polynomials satisfy a unique 
differential equation of the form 

oo 

ai{x)y^^\x) + xy"{x) + (a + 1 — x)y'{x) + ny{x) = 0, 

i=0 

where {ai(x)}^Q are continuous functions on the real line and {aj(x)}^^ are independent of 
the degree n. In [|| H. Bavinck found a new method to obtain the main result of |^]. This 
inversion method was found in a similar way as was done in |Q] in the case of generalizations 
of the Charlier polynomials. See also section 5 for more details. In 0] we used this inversion 
method to find all differential equations of the form 



MY,a^{x)y^'^{x)+NY,h,{x)y^'^{x) + 

2=0 2=0 
OO 

+ MiV^Q(x)y(')(x) + xy"{x) + (a + 1 -x)y'{x) + ny{x) = 0, (1) 

2 = 

where the coefficients {02(2;)}^;^, {bi{x)}'^^ and {cj(x)}^-^ are independent of n and the 
coefficients ao(x), 60 (x) and co(x) are independent of x, satisfied by the Sobolev-type Laguerre 
polynomials \ -L^'*^'^(x) > which are orthogonal with respect to the inner product 

L J n=0 

1 ^00 

<f,9> = TT7—^ I x''e-^f{x)g{x)dx + Mf{0)giO) + Nf{0)g\0), 

00 

n=0 



r(a + 1) Jo 



where a > —1, M > and > 0. These Sobolev-type Laguerre polynomials |L^''^'^(a;)| 
are generalizations of the generalized Laguerre polynomials |-L°'*^(j;)| . In fact we have 

L^'^'O(x) = L^^^'Hx) and L-'O(x) = L(f)(x). 



1 



In this report we will prove an inversion formula involving the classical Jacobi polynomials 
which can be used to find differential equations of the form 



i=0 i=0 i=0 

+ (1 - x^)y"{x) + [P-a-{a + P + 2)x] y\x) + n(n + a + /3 + l)y{x) = 0, (2) 



|16|. 



where the coefficients {aj(x)}^;^, {hi{x)}'^i and {ci{x)}'^i are independent of n and the 
coefficients ao(x), 6o(x) and cq{x) are independent of x, satisfied by the generalized Jacobi 
polynomials \P^'^'^^'^{x)\ . As an example we apply the special case /? = a of this 
inversion formula to solve the systems of equations obtained in 

The inversion formula for the Charlier polynomials obtained in [Q] (see also section 5) 
was also used in to find difference operators with Sobolev-type Charlier polynomials as 
eigenfunctions. 

In (I H. Bavinck and H. van Haeringen used similar inversion formulas to find difference 
equations for generalized Meixner polynomials. 

2 The classical Laguerre and Jacobi polynomials 

In this section we list the definitions and some properties of the classical Laguerre and Jacobi 
polynomials which we will use in this report. For details the reader is referred to Q, |12] and 

Lfi (x) > can be defined by 

J n=0 

n / j_ \ k 

L'nH-) = i:yi)'[l^_iyj^, n = 0,l,2,... (3) 
for all a. Their generating function is given by 

(1 - t)-"-' exp (-^) = £ LtHxr (4) 

n=0 

and for all n € {0, 1,2,...} we have 

D^Li-Hx) = {-iyLi"^\x),z = 0,l,2,...,n, (5) 

where D = — denotes the differentiation operator. The Laguerre polynomials satisfy the 
ax 

linear second order differential equation 

xy"{x) + (a + 1 — x)y'{x) + ny{x) = 0. (6) 

It is well-known that 

^-gM)'(:!:)4°'w. -0.1.2 
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This formula can easily be proved by using definition and changing the order of summation 
as follows 

= Ef"^"Ki:Vl)^K'l=^,n = 0,l,2,.... 
The classical Jacobi polynomials |Pn"'^^(x)| can be defined by 

L J 71=0 

Pi-«W = ^ (n + a + Zj+Dtl^+^t+^lj^-t ,, = 0.1,2.... (8) 

= 2-"|:(:!:)("r)(--i)v+ir-'.«=o.i.2.... 

for all a and /?. For all n G {0, 1, 2, . . .} we have 



D^P^-^^\x)='-^^±^^±^PtV''^Hxl^ = ii^^2,...^n. (11) 

The Jacobi polynomials satisfy the symmetry formula 

f("'^)(-:E) = (-irPf'-)(x), n = 0,1,2,... (12) 
and the linear second order differential equation 

(1 - x^)y"{x) + [/3-a-(a + /3 + 2)x] y'{x) + n(n + a + /3 + l)y{x) = 0. (13) 
Further we have for a + /3 + 1 > (compare with |]14|, page 277, formula (30)) 



2 y ^0 (a + /3 + fc + l)„+i V ^ ; 

This formula is much less known than formula ^ for the Laguerre polynomials. However, 
the proof is quite similar. In section 6 we will prove a much more general formula. 
We remark that (|l^ can be written in a more general form as 

A (-n)fc(a + (3 + l)fc(a + k + l)„_fc(a + f3 + 2k + l) («,/3) , 



fc=o r(a + /? + n + A: + 2) ^ ^""^ 

1 /l _ rr\ " 

, n = 0,1,2,..., (15) 



T{a + I3+1)\ 2 
which is valid for all a and (3. 
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Pn ' {x)\ and the Laguerre polynomials { Ln (x) \ are 

J n=0 I. J n=0 

related by the limit 

L^^\x) = lim P^°'^) fl - 5V (16) 



2x 

We remark that if we replace x by 1 3" (0)' multiply by /3" and let (3 tend to infinity 

in the complex plane along the halfline where a + /? is real and a + /3 + l>0we obtain (j^) 
by using (|l6|) and the fact that we have for all n G {0, 1, 2, . . .} 

+ k + l)n-k = {-if {^^ ^ n!, A; = 0, 1, 2, . . . , n. 

Finally we derive some formulas involving Jacobi polynomials which we will need in this 
report. By using definition (^) we easily obtain 

nP^/^)(x)-(n + a)Pi^f+^)(x) 

" {n + a + P+l)k{a + k + l)n-k / x-l X^ 

""h ^- ("-^)' V 2 ; + 

, , {n + a + f3 + l)k{a + k + l)n^k~i (x-l^^ 

-^" + "^?o (-^-1)' 

^ {n + a + f3 + l)k{a + k + l)n-k , , (x-\\^ 

E- 

fe=0 



Hence 

(x - l)Z)Pi"'^)(x) = nP^^)(x) - (n + a)Pi^f (x), n = 1, 2, 3, . . . . (17) 



By using the symmetry property ( jT^ ) we also have 

(x + l)Z)P^"''3)(x) =nPi"'^)(x) + (n + /?)Pi^t''''^(^), n = 1,2,3,.... (18) 
Another easy consequence of definition (|8|) is 

P("+^'^)(X) - P^^+i)(x) = PtX'''^'\x), n = 1,2,3,.... (19) 

If /3 = a we can combine the relations (p^, (|18|) and ( p!^ ) to find that 

2xP'Pi"'")(x) = 2nPi"'")(x) + (n + a)Pi°+^'"+^^(x), n = 2,3,4,... . (20) 

This formula turns out to be very useful in section 7. In that section we will also need the 
following well-known relation for ultraspherical polynomials (see for instance ]l6t) 

(n 



+ 2a + l)(n + 2a + 2)P^"+i'"+i) (x) -{n + a){n + a + l)P^^+''"+'^ (x) 

= 2(n + a + l)(2n + 2a + l)Pi'''")(x), n = 2,3,4,.... (21) 
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This formula can be proved by straightforward calculations as follows. As before we find by 
using definition (|8|) 

(a+i,a+i) p("+i,a+i) , 



(n + 2a + l)(n + 2a + 2)P^"+^'"+^^(x) - (n + a)(n + a + (x^ 

/ ^9 ^9 ^ (n + 2a + 3)fc (a + fe + 2)„_fc ( x - , 
(n + 2a + l)(n + 2a + 2)g (^^j + 

. , w , , ^,^^(n + 2a + l)fc(a + A; + 2)„_fc_2 /a;-!^^ 
-(n + a)(n + a + l)g (n-fe-2)! 

" (n + 2a + l)fc (a + fc + 2)„_fc 

X [(n + 2a + A; + l)(n + 2a + A; + 2) - (n - A:)(n - A: - 1)] 

9^9 ^9 , (n + 2a + l)fc(a + fc + l)n- fc+i ( x - ^ 

= 2(2n + 2a + l)> -. — - 

k\ {n-k)\ 

= 2(n + a + l)(2n + 2a + 1)^^"'"^ (x), n = 2, 3, 4, . . . . 
Another useful relation involving ultraspherical polynomials is 

(a + l)Pi°+i'°+i)(x) - (n + a + = i(n + a + 1)(1 - x^)^^^^ (22) 



which holds for n = 2, 3, 4, . . .. This can easily be shown by using definition ( [10| ) as follows 

(a + l)Pi"+i'"+i)(x) - (n + a + 
(« + l)l 



fc=0 L 



[X) 

'n + a + l\ /n + a + l' 
, n - k \ k 



-{n + a + l) 



n + a\ n + a 
n — k ) \ k 



n-l / 

-(n + a + l)2-" V 

— ^ \ 

fc=i \ 



n-2 



(x-l)'=(x + lf-'= 

(x-i)'=(x+ir-^ 

n + a W?^ + aV , -,^n-k-2 



( n + a \ 


( n 




\n-k-l] 


[k 





k=0 



= i(n + a + l)(l-x2)f'(^f'""''^(^), " = 2,3,4,.... 

3 Some summation formulas 

In this section we will derive some summation formulas which will be used in this report. We 
start with a partial sum of a iFq hypergeometric series at the point 1. In fact we have 

^(a + k\ " (o + l)^ (g + 2)„ (n + a + l\ 
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This formula appeared as a special case of lemma 4 in Q. Further we have the well-known 
Vandermonde summation formula 



2^1 



-n, b 
c 



(C - b)r 
(c)n 



(c)„ ^0, n = 0,1,2,.... 



(24) 



By using (|24[) we find that 



k=0 
6 2^1 



(6 + n + l)fcA;! 
/ -n,6 + 1 
\b + n + l 



nb 



■ 2F1 



b + n+1 
nb 



(6 + n + l)„ b + n + 1 {b + n + 2)n-i 
We remark that this can be written in a more general form as 

{-n)k{b)k 



-n + 1,6 + 1 
b + n + 2 

0, (6 + n + l)„/0, n = l,2,3,.... (25) 



^^T{b + n + k + l)k\ 



{b + 2k)=0, n = 1,2,3, ... , 



(26) 



which is valid for all b. This formula will be used in section 6. 
We also have the well-known Saalschiitz summation formula 



3F2 



-n,a,n + b + c — 1 
b,a + c 
Note that this implies that 



(6 - a)n(c)n 
{b)n{a + c)n 



, {b)n{a + c)n ^ 0, n = 0,l,2,.... (27) 



(c)n 



3-^2 



-n,^{b — c),n + b + c — 1 
thi 

-n)kib)k{c)k 



1 , (6)„((6 + c)/2)„ / 0, n = 0, 1, 2, . . . . (28) 



This can be used to show that for {^b + ^)n{b + n + l)n{b — c + 1)„ 7^ 



^,{b + n + l)k{b-c+l)t 



kl 



' + 2k) 



(6)„,+i(i6- c + i)„ 



c+l)„(i6 + i)„ 



(29) 



So we suppose that n e {0, 1, 2, . . .} and that {^b + ^)n{b + n + 1)^(6 — c + 1)„ / 0. Then we 
use (28) and change the order of summation to obtain 

{-n)kib)k{c)k 



E 



■'^{b + n + l)k{b-c+l)kkl 



' + 2k) 



i-n)k{b)k , p /-fc, i6-c + i,6 + A; 
(6 + n + l),fc!(^ + '^)^^^l, 6-c+l,l6+i 



E 

fc=0 
n 

^J6 + n + l)fcfc!' (6 _c+ 1)^(15 + i)^m! 



-n)k{b)h 



jib + 2k) J2 



{-k)rn{^b - C + l)„,{b + k)r 



n n~m 



{-n)m+ki-m - k)m{b)2m+k{\b - C + \)m 

tok^o + " + '^)m+k{b -c + l)m{^b + l)m{m + k)] ml 

_ i-n)mib)2mi ^b - C + 

m=0 



E E 

m=0 fc=0 

n 

E(-ir 



(6 + 2m + 2A;) 



{b + n + l)^{b-c+l)^{\b+\ 



2)m 



A;=0 



+ n + m + l)fc/c! 
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Hence, by using ( |25| ) we conclude that 

{-n)k{b)k{c)k 



(-ir 



-n)n{b)2n{\b 



c + 



(6 + n + l)„(6-c + l)„(i6 + i)„ 



-{b + 2n) 



{b)2n+l{\b- C+ \)r 



_ {b)n+l{\b - C+\)n 

{b + n+l)n{b-C+l)n{\b+\)n ~ {b - C + l)n{\b + \)n 

which proves (pT 



n = 0,1,2, 



4 The systems of equations 

Let a > — 1. The Sobolev-type Laguerre polynomials can be written as 

L J n=0 

L^'A^,7V(x) = AoL'^^\x) + ^iZ)4°)(x) + A2D^L'C\x), n = 0, 1, 2, . . . , 
where the coefficients and ^2 are given by 



n-iy (a + l)(a + 3) - 2 

MiV /n + a\ /n + a + 1 
n - 2 



+ 



(a + l)(a + 2) \n- ly 
n-l,yn + a\ 2MiV ( n + oi\ ( n + a + 



In ; a + l \n-l {a + lf 



n 



n-2 



N (n + a\ MN n + a\ n + a + l 



a + l\n — 1 (a + l) 



n 



n — 1 



For details concerning these Sobolev-type Laguerre polynomials and their definition the reader 
is referred to Q and [|ll|. Since the classical Laguerre polynomials |Ln"^(x)| ^ satisfy the 
differential equation (^) it is quite reasonable to look for differential equations of the form 
(|l]) for these Sobolev-type Laguerre polynomials |l"'^'^(2;)| ^ in view of this definition 

and the fact that L'^'^'^{x) = Li^\x). In it is shown that this leads to eight systems of 
equations for the coefficients {ai{x)}'^Q, {bi{x)}'^Q and {ci(x)}^Q. In order to find these 
coefficients we have to solve systems of equations which are of the form 



00 



i=l 



Ai{x)D'+''Ll^^x) = Fn{x), n = k + l,k + 2,k + 3,..., 



where /c € {0, 1, 2, . . .} and the coefficients {Ai{x)}°^^ are independent of n. In |^ it is pointed 
out that this system of equations has a unique solution given by 



A,{x) = i-iy+>'J2L\-_; 



(-a-i-fc-l) 



{-x)Fj+kix), i = 1,2,3,.... 
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oo 





This is an easy consequence of the Laguerre inversion formula 

jZLtV''\-^)Lt^p{x)=5i,, j<i, i,j = 0,l,2,..., (30) 

k=j 

which was found by H. Bavinck in For more details the reader is referred to and Q. 
See also section 5 of this report. 

Now we take a > —1 and /3 > — 1. The generalized Jacobi polynomials 
can be written as 

pa,/3,M,N^^^ = AoPt'^^Hx) + [Ai{l - x) - ^2(1 + x)] DP^^^'^^x), n = 0, 1,2, . . . , 
where the coefficients and A2 are given by 

fn + p\ fn + a + P + l\ fn + a\ fn + a + P + 1\ 

In— 1/1 n I In— 1/1 n / 
Ao = l + ^ ^ + ^ ^ + 

+ M7V-^" + ^ + ^^' /^n + a + /3 + lV 



(a + l)(/3 + l)l n-1 



n + /3\ / n + a + /3^ 



^_ M ynyy n y MiV /n + a + /3\/n + a + /3 + l 
^ a + P + l /n + a\ a + ly n - 1 J y n ^ 



(n + a\ (n + a + I3\ 

N \ n j\ n j MN fn + a + p\fn + a + P + l' 
Ai = : 7 r h 



a + /? + 1 









/3 + l\ n — 1 I \ n 



Here we used the same definition as in []13| , but in a slightly different notation. The case 
a + P + 1 = must be understood by continuity. In view of this definition and the fact 
that the classical Jacobi polynomials satisfy the differential equation (^) it 

is quite natural to look for differential equations of the form (|2|) satisfied by these generalized 
Jacobi polynomials as was already pointed out in Again this leads to eight systems of 
equations for the coefficients {ai(x)}^Q, {bi{x)}^Q and {cj(x)}^Q. In order to find these 
coefficients we have to solve systems of equations which are of the form 



A{x)D'Pt'''Hx) = n = 1, 2, 3, . . . , 



4 = 1 

I 00 



where the coefficients {Ai{x)}°^i are independent of n. This system of equations has a unique 
solution given by 



MX) -2 [x)I'j{x),i-L,2,6,.... 
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This is a consequence of the Jacobi inversion formula 



E 

k=j 



a + (3 + 2k + l 



ia + P + k + j + 



i—k 



{a+j,f3+j) 
k-j 



(x) = 6ij, j <i, i,j = 0,1,2,... , (31) 



which will be proved in this report. Again, the case a + /3 + 1 = must be understood by 

2x 

continuity. We remark that if we replace x by 1 in (31), multiply by (3^~^ and let (3 tend 

/5 ^ 



to infinity along the positive real axis we obtain the Laguerre inversion formula (^0|) by using 

(H). 

In pH] we found all differential equations of the form 



M^ai(x)y(*)(x) + (1 - x^)y"{x) - 2{a + l)xy'{x) +n{n + 2a + l)y{x) 

i=0 



(32) 



where {aj(x)}^Q are continuous functions on the real line and {ai(x)}^]^ are independent 
of n, satisfied by the symmetric generalized ultraspherical polynomials P"'"'^^'^^(x) ^ 

L J n=0 

defined by 

pa,a,M,M(^x^ = CoPt'^K^) " CixDPt'''\x), n = 0,1,2,..., 



where 



2Mn n + 2a + l\ _,n/n + 2Q + l' 

1 H + 

a + 1\ n / \n— 1, 



2M (n + 2a\ 2K-P /n + 2a\ /n + 2a + 1^ 
2a + 1 \ n / a + 1 \ n-1 ]\ n 



We remark that these polynomials form a special case {(3 = a and N = M) of the generalized 
Jacobi polynomials |p°'^'^^'^(x)| ^, but the differential equation (|3^) has a very special 

form without a M^-part. This differential equation will appear not to be a special case of the 
differential equation of the form (Q) for the generalized Jacobi polynomials, since the MN- 
part will not vanish if we take (3 = a and = M. We aim to give a proof of this in a future 
publication. In section 7 of this report we will apply the special case /3 = a of the Jacobi 



inversion formula (|3l| ) to solve the systems of equations obtained in |1C] as an example. 



5 The inversion formulas 

In Q H. Bavinck and R. Koekoek found the following inversion formula involving Charlier 
polynomials 



E Ct^i-x)d^}^{x) = 5ij, j <i,i,j = 0,1,2,.... 

k=j 

This formula is an easy consequence of the generating function (see for instance |T2| 

oo 

e-^^l + tf = ^C7^")(x)r. 

n=0 



(33) 
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In fact we have 



k=0 m=0 



CO I n 



n=0 \fc=0 / 



Hence 



k=0 



1, n = 



0, n = l,2,3,.... 



Now (33) easily follows by taking n = i — j and shifting the summation index. This formula 
was also used in to find difference operators with Sobolev-type Charlier polynomials as 
eigenfunctions. In Q a similar formula involving Meixner polynomials was used to find 
difference equations for generalized Meixner polynomials. 

Formula (33) can be interpreted as follows. If we define the matrix T = {tij)'^j^Q with 
entries 



ti 



.0, j> i, 

then this matrix T is a triangular matrix with determinant 1 and the inverse U of this matrix 
is given by = U = {uij)i'j=o with entries 

Ctf{-x), j<i 



0, 



j > i. 



Therefore we call ( p3D an inversion formula. 

In the same way we find by using the generating function (Q) for the Laguerre polynomials 



k=j 



i-a-2). 



Sij, j < i, i,j = 0,1,2,, 



However, this formula cannot be used to solve systems of equations of the form 

oo 

^ ^,(x)D^L(f ) (x) = Fn{x), n = 1,2,3,... 



i=l 



in view of the parametershift in (|5|) . 

In [|l| H. Bavinck used a slightly different method to find the Laguerre inversion formula 
(30) from the generating function (^) for the Laguerre polynomials. In fact we have 



(1-t) 



(l-t) 



~a-j-l 



exp 

oo oo 



xt 

t - 1 



(1 - exp 



—xt 

t - 1 



/c=0 m=0 



oo / n 



E E4"^^^(-)4i. 

n=0 \A:=0 
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This implies, by comparing the coefficients of f ^ on both sides, that 



k=0 



which is equivalent to (|30|). 

Formula (30) can be interpreted as follows. If we define the matrix T = {tij)'^j^Q with 
entries 



t 



0, j > i, 

then this matrix T is a triangular matrix with determinant 1 and the inverse U of this matrix 



IS given 



by T~^ = U = (u,,) 



i3Ji,j=Q 



with entries 



0, j > i. 

It is possible to generalize the Laguerre inversion formula to 



k=0 



{Pn - qn + 2)n 



n! 



, n = 0,l,2, 



(34) 



where Pn and qn are arbitrary and even may depend on n. In order to have an inversion 
formula we have to choose p„ and g„ such that 

{Pn - g„ + 2)„ = 0, n = 1, 2, 3, . . . , 

hence 

Pn-Qn^ {-n - 1, -n, . . . , -3, -2}, n = 1, 2, 3, . . . . 

Note that the endpoint-cases Pn — Qn = —2 and Pn — In = —n — 1 correspond to the earlier 
mentioned inversion formulas. 

To prove (U) we use (§) to obtain 



(1-t) 



-a-l 



exp 



xt 

t - 1 



n = 0,l,2,... 



where Dt = — denotes differentiation with respect to t. Hence by using Leibniz' rule we find 
dt 

n 



ELr-^(x)LrT^"^(--) 

fc=0 



n , 



fc=0 



(1-t) 



-a-pn-l 



exp 



xt 
t-l 



t=0 



-A" (l-i) 
n! L 



(n - A:)! 

qn-pn — 2 



n~k 



\a+qn-l , 



(1 - ty ' ^ exp 

{Pn -qn + 2)„ 



-xt 



t - 1 

, n = 0,l,2. 



i=0 



i=0 n! 
which proves (34). 

In case of the Jacobi polynomials the above methods seem not to be applicable. However, 
in the next section we will give a proof of the Jacobi inversion formula (^) . 
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6 Proof of the Jacobi inversion formula 



In this section we will prove that 



^ {a + P + 2k + l){a + P + l)k p(a,/3)^ /-n-a-l-n-P-l) , . 

to na+p+n+k+2) ^''^^"-'^ 



1 1 f X — y 



r(a + /3 + l)n! V 2 



n = 0,1,2,..., (35) 



which holds for all a and /?. 

Note that (|l5|) is a special case of (|35|) since 



(_n_a-l-n-/3-l). . ^ (-n-a)n_fc ^ .„_fc (Q + fe + l)n_fc 

^ ^ (n-A:)! ^ ^ {n - k)l 

(-ir, 



n! 



-(-n)fc(a + A; + l)„_fc, A; = 0, 1, 2, . . . , n 



for all n G {0,1,2,...}. 

By taking y = a; in (^5|) we easily obtain 



^ (a + /3 + 2fc + l)(a + /3 + l)fc 



T{a + p + n + k + 2) 

i n = 

r(a + /? + l)' ^3g^ 



0, n = l,2,3, 



for all a and p. If we take n = i — j in (36) and shift the summation index we find 

A {a + P + 2k-2j + l){a + P + l)k-j ^ 
t T{a + p + i-2j + k + 2) ^ 

K—J 



For a and /3 real with a + /3 + l > —1 we now obtain ( |3lD by shifting both a and /3 by j. 
Note that (]35| ) for y = — x in a similar way leads to 



^ a + p + 2k + l 



(a + /3 + A; + j + 1) 



X '^'■^(^) = J- < ^> = 0, 1,2, ... , (37) 



which will also be used in the next section. 

In order to prove ( ^5| ) we start with the left-hand side, apply definition (|8|) to P^'^\x) 
and definition (^) to a i, n p i) change the order of summation to obtain 

^ {a + P + 2k + l)(a + P+l)k p{a,l3) , .p{-n-a-l-n-(3-l), . 

t T{a + P + n + k + 2) ^ ^""l^n-k ^V) 
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" ^ (g + /3 + 2fc + + /3 + l)fc (g + /? + fc + l)i (g + i + 1 ) 

T{a + (3 + n + k + 2) i\ {k - i) 



A;=0 1=0 j=0 



k—i 



{a + P + n + k- j + 2)j {a + k + l)n-k~j f x - IV ( y - IV 



j\ {n-k-j)\ 

n n n—k 
1=0 k=i j=0 

{a + [3 + 2k + l){a + P + l)i+k{a + i + l)n-t-j f x - IV fy-lV 



T{a + P + n + i + k - j + 2)i\k\j\ {n - i - k - j)] \ 2 

n n—i n—i—j , i\i 

^ X — \\ I y — ly 



r(a + /3 + n + /c - j + 2) i! (/c - i)\ j\ {n - k - j)l \ 2 

n n—i n—i—k 

EE E (-ir-^-^x 

i=Ofc=0 j=0 

(g + /? + 2i + 2fc + l)(a + /3 + l)2^+k{a + i + l)n-i-i ( x - l\ (y-lV 



EE E (-ir 

i=Oj=0 fc=0 



X 



(g + /? + 2i + 2A; + l)(a + /3 + l)2*+fc(g + i + l)„-»-j(-n + i + j) 



r(a + /3 + n + z - j + /c + 2)i! j! A;! (n - i - j)! 



•^Y^f ,„_j (g + /3 + l)2»(g + i + l)«-i-j / a: - l y - 1 
^^o.t^o i!j!(n-.-i)! I 2 M 2 



X 



(-n + i + j)fe(a + /3 + 2i + l)fc ^ 
X 2^ . ^ . — —. — T—iy^ g + P + 2z + 2A; + l , n = 0,1,2,.... 

Now we use ( p6D to obtain 

^ {a + P + 2k + l)(a + /? + 1)^ ^{a,l3) ^ x pC-n-a-l . . 

^t^Q r(a + /? + n + A; + 2) ^"^'^^"^ 
^(_l)"-(" + '^ + l)2i ^^-l^Vy-l^""' a + /? + 2i + l 



i=0 



il{n-i)\ \ 2 J \ 2 J T{a + p + 2i + 2) 



1 l|,/nU-lV/-l-. 



r(a + ^ + l)n!^^ V 2 / V 2 

n = 0,1,2,..., 



r(a + /3 + l)n! V 2 
which proves (|35|). 



7 The symmetric generalized ultraspherical polynomials 



Let g > — 1. In [T^] we found the coefficients {aj(x)}^Q of the differential equation ( p2| ) for 
the symmetric generahzed ultraspherical polynomials \ P"'"'^^'^^(x) > . In order to do this 

L J n=0 
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we had to solve the following two systems of equations for the coefficients {ai{x)}^ 



i=o ■ 



and 

oo oo / I Q I 1 \ 



i=0 1=0 



for n = 0, 1, 2, . . ., where the coefficients {aj(x)}^Q are continuous functions on the real line 
and {0^(2;)}^^ are independent of n. Now we suppose that 00(2;) := ao{n,a) is independent 
of X as we did in 0. Then it is clear (see for instance lemma 1 in 0) that ai{x) must be 
a polynomial in x of degree at most i for each i = 1,2,3, .. .. In |jl^ we showed that the 
solution for {0'i{x)}^Q is not unique. In fact it was shown that 

ao{x) := ao{n, a) = ao(l; Ci)bQ{n, a) + co(n, a), n = 0, 1, 2, . . . (40) 

and that 

ai{x) = 00(1,0)64(2;) + Ci{x), f = 1,2,3, ... , (41) 
where ao(l,a) is arbitrary and 

5o(n, a) = ^[l- i-ir] , n = 0, 1, 2, . . . , (42) 
co(n,a) = 4(2a + 3)( j , n = 0, 1, 2, . . . , (43) 

h{x) = '^{-xy, i = 1,2,3,..., (44) 

ci(x) = and Ci{x) = {2a + 3)(1 - x^)'^P^^^-'^^'''-'^^\x), i = 2,3,4, ... . (45) 

In this section we will give an alternative proof of this by using the Jacobi inversion 
formula (|3l|). Note that the special case /3 = a of the Jacobi inversion formula (|3l|) reads 

2a + 2k + l 
f^.j2^TkTJTT)~J7i 

X /^^-r^-^'— '-'^(x)p£+^'"+^)(2;) = J < z, ^,J = 0, 1,2, ... . (46) 

If we apply this inversion formula to the system of equations 

00 

J2Mx)D'Pt''\x) = Fnix), n = 1,2,3,..., (47) 

i=l 

then we find 

Mx) = P/_-— ^^(x)F,(x), i = 1,2,3, ... . (48) 
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The special case /? = a of (^) reads 

^ 2a + 2fc + 1 
^. (2a + fc + j + l),_,+i 



X P. 



(— a— j— 1,— a— i— 1) 



i—k 



, j<i, i,j = 0,1,2,.... (49) 



By considering ( [3^ ) and (^) for n = and n = 1 we conclude that ao(0, a) = 0, ao{l,a) 
is arbitrary and 01(2;) = —ao{l,a)x. For n = 2, 3, 4, . . . it turns out to be more convenient to 
use another system of equations instead of (^). By using (^) we find for z = 0, 1, 2, . . . 

= nD^Pt''\x) + l{n + a)D^Ptf'''^'\x), n = 2,3,4,.... 



Combining (^) and (|39|) we now obtain 

00 r) 



n + a 



n + a 

8 / n + 2q 



/n + 2a + l\ 4n /n + 2a\ 
y n - 1 y ~ 2a + 1 \ n J 

n + 2a + l\ (n + 2a 
n — 1 I In— 1 



n + a\ n — 2 / 

So we conclude that ( |39|) for n = 2, 3, 4, . . . may be replaced by 



x), n = 2, 3, 4, . . . . 



ip(a+l,a+l), 



X) 



i=0 



n + a\n-2 j " 



n = 2,3,4,.... 



(50) 



Note that for n = 2 this implies that ao(2, a) = 4(2a + 3). 

Since ai{x) must be a polynomial in x of degree at most i for each i = 1, 2, 3, ... we may 
write 

ai[x) = kiX^ + lower order terms , i = 1, 2, 3, . . . . 
By comparing the coefficients of highest degree in ( |3^ ) and (^) we find by using (|8|) : 



ao(n,a) _^ ^ fcj 



ni 



(n — iV. 



0, n= 1,2,3,... 



and 



ao(n, a) 



n-2 



+ E 



4(2n + 2a - 1) 



n + 2a\ 1 



n = 3, 4, 5, , 



(n-2)! ^(n-i-2)! ' ' \n -2 j (n 

Since ki is independent of n for i = 1, 2, 3, . . . and ao(2, a) = 4(2a + 3) we conclude that 

ao(n + 2, a) - ao(n, a) = 4(2n + 2a + 3) , n = 0, 1, 2, . . . , (51) 
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where 00(0,0;) = and ao(l,a) is arbitrary. Hence 



ao(2n,a) -ao(0,a) = 1^ ) (4A; + 2a + 3), n = 1, 2, 3, , 



and 

n-l 

2k + 1 



oo(2n + 1, a) - ao(l, a) = 4 ^ ( " + ^ ) {4k + 2a + 5), n = 1, 2, 3, . . . . 



fe=0 

Since we have by using (p^) 

2n-2 



fc=0 

and 

n-l 



^ (2a + 3)fc+i , ,, (2a + 5)2n-2 /2n + 2a + 2\ 

g = + (2n-2)! =(^" + 3)(^ 2n-2 j'- = l'2'3, 



^,2/c + 2a + 3\ (2a + 3)2^+2 , ^ (2a + 3)2^+1 

2^^(2a + 3)fc+i , ,, (2a + 5)2n-i /2n + 2a + 3\ 

£ = + (2n-l)! =(^« + 3)(^ 2n-l j'- = l'2'3, 



we now conclude that (|4g), (||) and (gSj) hold. 
The systems of equations (|3^ ) and (^) lead to 



£a,(x)Z)^P^")(x) = _^(^^+^2"jz)2pK")(x) -ao(n,a)P^")(x) (52) 
for n = 0, 1, 2, . . . and 

f a, (x)D^Pi"_t^'"-^^) (x) = -A- f " Z)2p("'") (x) - ao (n, a)PS''"^'^ (x) (53) 

1=1 \ / 

for n = 2, 3, 4, . . .. 

First we remark that (52) is true for n = and n = 1 since ao(0, a) = and ai(x) = 
— ao(l,a)x. Then we will show that every solution of ( p^) also satisfies ([5^). Suppose that 
{oi(x)}^]^ is a solution of (^). Now we use (0), (^T]), (^iD and the fact that {ai{x)}°l^ are 
independent of n to obtain for n = 2, 3, 4, . . . 

00 

2(n + a + l)(2n + 2a + 1) ^ a,(x)P»*P^"'") (x) 

i=l 

00 

= (n + 2a+ l)(n + 2a + 2)^ai(x)D^Pi"+^'"+i)(x) + 

00 

- (n + a) (n + a + 1) ^ a, (x)Z^'P^"+''"+') (x) 
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(n + 2a + l)(n + 2a + 2) 



8^ /n + 2a + 2\ , , _ ^^^^ ^ ^ 

n + a + 2 \ n 



+ 



(n + Q)(n + a + 1) 



8 fn + 2a^ 



n + a \ n — 2 
2{n + 2a + l)(n + 2a + 2) ( n + 2a + 2\ 



D^Pt^-){x)-a,{n,a)PtX'''^^'\x) 



n + a + 2 



n 



-8(n + a + l) 



n + 2a' 
, n-2 , 



+ 2a + 3)(n + 2a + 4)Pi"+2'"+2) (x) + 

-2(n + a + 2)(2n + 2a + 3)P^"+^'"+^) (x^ 

- ao(n, a) [(n + 2a + l)(n + 2a + 2)P^"+i'"+^)(a;) + 

-(n + a)(n + a + l)Pi^+''"+') (x)" 

2(n + a + l)(n + 2a + l)(n + 2a + 2) (^"^ + + P£t''"^'H^) + 



+ 



(n + a + 1) 



n + 2a^ 
n-2 , 



2(n + a + l)(2n + 2a + l)ao(n,a)P^"'")(x 
n + 2a + 2\ fn + 2a 
n J \ n — 2 

2(n + a + l)(2n + 2a + l)ao(n, a)P^"'") (, 

4 /n + 2a' 
2a + 1 I n 



(n + a+ 1) 



2(n + a + l)(2n + 2a + l) 



P'2p^°-°)(x) - ao(n, a)P^°'")(x) 



Since a > — 1 this proves that every solution of ( |53D also satisfies 

Now we will solve (^). Shifting n by two we may write, since the coefficients {aj(x)}^;^ 
are independent of n 



X: a,(x)Z)*Pi"+i'-+i)(x) = F„(x), n = 0, 1, 2, . . . , 



(54) 



1=1 



where 



Fn{x) 



+ + ^ j Z)2pi:'? (x) - ao(n + 2, a)p(-+i'"+i) (x) . 



n + a + 2 

Since ao(2,a) = 4(2a + 3) we easily find that ^0(2^) = 0- This implies that the system of 



equations ( |54|) is of the form (47). So if we apply the inversion formula ([4q) to the system of 
equations (|53) we obtain by using (p 



a,: 



(^■^ _ oi V 2a + 2j + 3 

^ ^"\ei(2« + J+3)m^-^' 



(x)F,(x), i = 1,2,3,.... 
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Hence, by using ( |40| ) we now conclude that the coefficients {aj(x)}^| can be written in the 
form (|l|). Moreover, we find by using (H), (|l2|) and (||) 

Hx) = 2-^ ± (^IV/I t)L ^17""''°'^"^ (x)P,^°+^'"^^) (X) [(-1)^- - 1 

_ ni-l 2a + 2j + 3 (_„_i_2 -a-i-2) . , p{a+l,a+l) , 

,eo(2« + J+3)m^-^' ' ^ 

^tj2a+j+3).+/-^- 

_ pi— 1 (~-^) • 1 r) q 

— z , z — i, z, o, . . . , 



which proves (Q). And by using (^3|) and ( [Tl| ) we obtain 

Ci{x) = 2 2^ i^Lj ^(x)Gj(x), ? = 1,2,3,..., 

where 

4(2a + 3) fj + 2a + 4\ 



G,(x) 



j+a+2\ j 



(a + 2)pf (x) - (j + a + 2)P,("+^'"+^)(x) 



It is clear that Gi{x) = 0, which implies that ci(x) = 0. Note that since = —x this also 
implies that ai(x) = — ao(l, a)x, which agrees with what we have found before. Now we use 
dH) to find 

G,{x) = (2a + 3)(1 - x') + + -^^^Cx), j = 2, 3, 4, . . . . 

Hence, for i = 2, 3, 4, . . . we have 

Ci{x) = (2a + 3)(l -x2)2^ x 

^ 2a + 2j + 3 /j + 2a + 4\ (-a-i-2,~a-i-2) , . p(a+3,a+3) . ^ 

"^^j2a + , + 3Wil, J j^-^- ^"^^^-^ 
Now it remains to show that 



^ 2a + 2j + 3 /j + 2a + 4\ fa-i~2-a~i~2). . p(a+3,a+3) 



_ 1 

In order to do this we write for i = 2, 3, 4 



_^(.-.+3,«-.+3)^^)^ ^^2,3,4,.... (55) 



^ 2a + 2j + 3 (j + 2a + A\ (_c.~i-2,-a-i-2) , . (a+3,a+3) . X 



i=2 

i-2 



E ^u: -r z.ft -r I ^^^Lt -r p(-a-^-2,-a-^-2) . x p(a+3,«+3) / n 

(2a + A: + 5),+i {k + 2)\ ^-^-^ ^'^^^'^ 
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Now we apply definition (^) to P^"^^'""*"^^ (x) and definition (^) to Pj^_^_2 ^' " * and 
cliange the order of summation to obtain for i = 2, 3, 4, . . . 

2a + 2fc + 7 (2a + 5)fc+2 p(~a-i-2,-a-i-2) . ^p{a+3,a+3), n 

^^j2a + fc + 5),+i (A: + 2)! ^-^^'^ ^"^^^^ ^''> 

i~2 i-k-2 k 
fc=0 m=0 n=0 

^ (2a + 2fc + 7)(2a + 5)fc+n+2(« + n + 4),„^_n-2 /x-lY'"+" 
^ (2a + A; + 5)i_^+i(fc + 2)!m!(f - A:-m-2)!n!(fc-n)! V 2 J 

i—2 i—m—2i—m—n—2 
m=0 n=0 fc=0 

^ (2a + 2n + 2fc + 7) (2a + 5)2n+fc+2(a + n + 4)^_^_„_2 /x^ix 
(2a + n + A; + 5)i_m+i(ri + A; + 2)!m! (i - m - n - A; - 2)!n! A;! V 2 

= EE E (-ir-^x 

j=On=0 fe=0 

^ (2a + 2n + 2fc + 7)(2a + 5)n+fc(a + n + 4)^_j_2 /x - ly 

^ (2a + 2n + A: + 7)i_j_i(n + A; + 2)! (j - n)! (z - j - A: - 2)!n! A;! V 2 

VVf-iy-" (2a + 5)„(a + n + 4),_,-_2 / x - ly 

,4^o^o (2a + 2n + 7)i_j_i(n + 2)!(j-n)!(i-j-2)!n! I 2 

Hence, by using (|29|), (|2^) and definition (^) we obtain for i = 2, 3, 4, . . . 

i-2 



2a + 2k + 7 (2a + 5)^+2 jy{~a-i~2-a~i~2) p(Q+3,a+3) ^ 

^J2a + + 5),+i (A: + 2)! 

'"^ ^' (2a + 5)„(a + n + 4)i_j_2 /x-lV' 



^ ^u^ -r z,A. -r I v^u^ -r ^;fc+2 (-a-»-2 -a-^-2) , x „(. 

ro., ^ i. ^ ^) .^^ (A; + 2V ^-^•"2 ^-^^^fe 



,t^o^o (2a + 2n + 7)i_,_i(n + 2)!(j-n)!(i-j-2)!n! V 2 

^ (2a + 2n + 7)^-^-1 (-a - l)i-j-2 
(n + 3)i_j_2(a + n + ^)i-j-2 

V ly-" (2a + 5)„(-a-l),,,„2 /x-iy 
t^.-n (n + i-j)!(j-n)!(i-i-2)!n! I 2 



X 



a — l)j_j_2 /^x — ly ^ j,2a + 5 



J=0 71=0 

i-2 

^ ^ ,eo(^-i)K^-i-2)!j! I 2 ; ' H + 1 

(-1)^ ^ (-g - l)i-j-2(i - j - 2a - i)j fx-l\^ 
i^- Po {i-j-2)\j\ 

(-1)'"^ ^ (^ - j - 2a - 4)j {-a - l)i-j-2 ( X - l y _ (a_i+3,a_i+3) 

^\ j^^ j\ {^-J-2)l [ 2 ) "zl^^^ 
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which proves (55). 



8 Some remarks 

Note that we have from definition (^) for the Laguerre polynomials that 

4-")(x) = (-ir^, n = 0,1,2,.... (56) 
n! 

Hence, the polynomial Ln ^\x) reduces to a monomial of degree n for all n € {0, 1, 2, . . .}. 
Definition ( p!o| ) for the Jacobi polynomials leads to 

which is also a monomial. However, this monomial might reduce to the zero polynomial. For 

instance, Pn ^\x) equals the zero polynomial for all n G {1, 2, 3, . . .}. 

2x 2y 

Further we remark that if we replace x by 1 — and y by 1 — in (35), multiply by 

P P 

r(a + /3 + l)/3" and let (3 tend to infinity in an appropriate way we obtain by using ( p^ ) 

E Lt\x)L^---'\-y) = n = 0, 1, 2, . . . . (57) 

fc=o ^■ 

Note that (0) is a special case of ( p7[ ) since 

4—- "(0) = (-i)-' '°|„'_\')'r'' = (:!;). * =0,1,2..... „ 

for all n G {0, 1, 2, . . .}. Moreover, note that ( p7[ ) is a special case of the well-known convolu- 
tion formula for the classical Laguerre polynomials 

n 

E Lt\x)L^^lM = Kr^^'Kx + y), n = 0, 1,2, . . . 

k=0 

in view of (|5^) . By using the technique demonstrated at the end of section 5 this convolution 
formula can be proved for all a and /? which might even depend on n. 
We remark that, by using the fact that 

{b/2Ub + 2k) = 6(5/2 + l)fc. A: = 0, 1, 2, ... , 



the formulas ( pq ) and ( p9| ) can also be obtained by using summation formulas for terminating 
well-poised hypergeometric series (see for instance formulas (HI. 9) and (III.26) in [p^]). 

The computation of oo(n,a) from ([5l] ) can also be done in the following way. Note that 
we have 



(2n + 2a + 3) ( " + + ^ j = (2a + 3) 



n + 2a + 4\ /n + 2a + 2' 
n / \ n - 2 , 



, n = 0,l,2,, 
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Hence, instead of ( pS] ) we can also use the telescoping property of the sums to find that 

'2n + 2a + 2^ 



and 



n-1 

E 

k=0 
n-1 
k=0 



2k + 2a + 2' 
[ 2k , 

'2k + 2a + 3\ 
2A; + 1 / 



{Ak + 2a + 3) = (2a + 3) 



2n-2 



{Ak + 2a + 5) = (2a + 3) 



2n + 2a + 3^ 
2n - 1 , 



I, n= 1,2,3, ... 



, n = 1,2,3, ... . 



Finally, we remark that if we apply the inversion formula ( ^6|) to the system of equations 
(52) instead of ( |5^ ) we find for i = 1, 2, 3, . . . that 



h.{x) = 2!-'Yl 



1 2a + 2j + 1 



(x)pf'")(x) (-1)^-1 



(58) 



and 



CiXx 



_L_J^- + 2")z..^K^)(,)_(2a + 3)| 



i + 2a + 2\ ^ 

J -2 ^""^ 



(59) 



From (]58| ) we easily obtain (^) in the same way as before by using (|l^), ( ^6|) and (49). 
Further we easily find from ( ^9[ ) that ci(x) = 0, but we were not able to derive ( p5[ ) for 
i = 2,3,4,... from (||). 
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